Entanglement Witnesses and Measures for Bright Squeezed Vacuum 



(N 

o 

(N 

> 
O 

(N 



> 
On 

o 

(N 



X 



Magdalena Stobinska, 1,2 Falk Toppel, 3,4 Pavel Sekatski, 5 and Maria V. Chekhova 3,6 

1 Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/4-6, 02-668 Warsaw, Poland 
2 Institute of Theoretical Physics and Astrophysics, University of Gdansk, 80-952 Gdansk, Poland 
3 Max Planck Institute for the Science of Light, Gunther-Scharowsky-Strafle 1/Bau 24-, 91058 Erlangen, Germany 
^University of Erlangen-Niirnberg, Staudtstrafie 7/B2, 91058 Erlangen, Germany 
Group of Applied Physics, University of Geneva, 
Chemin de Pinchat 22, CH-1211 Geneva, Switzerland 
6 Department of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory, 119991 Moscow, Russia 

Quantum entanglement is a fascinating phenomenon, especially if it is observed at the macroscopic 
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I. INTRODUCTION 

Quantum entanglement is a fascinating phenomenon, 
especially if it is observed at the macroscopic scale p], Q . 
It emerges from the quantum superposition principle ly- 
ing at the heart of quantum mechanics. If a two- or multi- 
mode superposition is highly populated, macroscopically 
entangled subsystems are created. A bipartite system 
maximally entangled in a given degree of freedom has 
two subsystems, for which the values of the degree of 
freedom are completely random but still perfectly cor- 
related. Importantly, macroscopic quantum correlations 
can be revealed only by accurate measurement outcomes 
and strategies Q. 

Entanglement is the basic resource for quantum infor- 
mation processing, quantum communication and other 
quantum technologies [3, [Ej]. F° r this reason, there is 
an important need for efficient and reliable entanglement 
verification and quantification in various physical sys- 
tems. This is done by the measurement of entanglement 
witnesses or measures ■ While an entanglement witness 
only tells whether a system is entangled or not, an en- 
tanglement measure allows one to quantify the amount 
of entanglement. Unfortunately, most of the entangle- 
ment measures are not operational, i.e. are not directly 
measurable in experiment. An operational measure of 
entanglement has been proposed for continuous variables 
such as wave vector or frequency Q, but it is absent for 
discrete variables of large dimension, e.g. 10 10 , such as 
photon number. 

The first condition for entanglement (inseparability), 
the Peres-Horodecki criterion, has been formulated in 
terms of density matrices || and therefore was not oper- 
ational. Later, other conditions were formulated in terms 



of measurable quantities such as variances or uncertain- 
ties of (9l-[l^| . A general method of finding inseparability 
criteria for bipartite continuous variable systems, deeply 
related to the Peres-Horodecki condition, was proposed 
in Q and then amended in (l3l . It contains a hi- 

erarchy of inequalities for the measurable moments of 
creation and annihilation operators and generalizes the 
previously obtained conditions j^-ll2j. 

Recently, macroscopic states of light generated a lot of 
interest and posed important questions in the scientific 
community. A possibility of verifying macroscopic entan- 
glement became intriguing and was widely discussed flB - 
Il7| . Among macroscopic photonic states, two kinds of 
entangled states have been generated and reported: the 
micro-macro polarization singlet state, where a macro- 
scopic qubit is entangled with a single photon [l8j . and 
entangled bright squeezed vacuum (BSV) states, macro- 
scopic analogs of two-photon polarization Bell states [19( . 
An entanglement test with the micro-macro singlet state 
was performed [18], but later refuted [20]. It was shown 
that in the case of macroscopic states, inefficient detec- 
tion may falsely reveal entanglement in separable states 
[2l| and that it is incapable of grasping their quantum 
character |22]. New criteria for the test were formu- 
lated [20L l23j|. but not put in experiment. Entanglement 
in BSV states is known as twin-beam multiphoton en- 
tanglement [13] and manifests itself in perfect polariza- 
tion correlations between two macroscopically populated 
beams. Macroscopic singlet polarization Bell state [l9| 
is formally equivalent to the singlet state of a large spin, 
with the spin value being hugely uncertain. An entangle- 
ment witness in the form of spin inequality for BSV was 
discussed |25l-[27||. Recently, the non-separability witness 
derived in Ref. [26| was tested in experiment [28| for this 
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state. The form in which the non-separability witness 
was derived is inapplicable to the other three BSV states 
(macroscopic triplet polarization Bell states) experimen- 
tally obtained in Ref. fl9] . Experimental realization of 
entanglement measures, e.g. the Schmidt number or other 
measures based on the eigenvalues of the density opera- 
tor, for entangled macroscopic states of light remains an 
open issue. 

In this paper, we formulate feasible non-separability 
(entanglement) witnesses for the macroscopic triplet po- 
larization Bell states in the spirit of the Duan criterion. 
Their testing involves efficient photodetection and the 
measurement of the Stokes operators variances. We also 
calculate measures of entanglement for these states such 
as the logarithmic negativity and the Schmidt number as 
well as the operational measure R introduced in Ref. [7|. 
Finally, we examine the dependence of the measure R and 
its relation to the Schmidt number on the experimental 
conditions. 

This paper is organized as follows. We start with in- 
troducing the basic properties of the bright squeezed vac- 
uum in Section|nl Section lHll is devoted to the derivation 
of the operational entanglement witnesses for the macro- 
scopic triplet polarization Bell states. Next, in Section llVl 
we discuss various entanglement measures: the Schmidt 
number, the effective Schmidt number, and the negativ- 
ity. A possibility of their experimental verification is ad- 
dressed. Finally, we make the conclusions in Section IVl 



II. ENTANGLED SQUEEZED VACUUM 

Entangled four-mode BSV states of light are pre- 
pared in experiment by employing two two-mode BSV 
states [1^1 obtained via high-gain unseeded parametric 
down-conversion (PDC). In Ref. they were created 
by overlapping on a polarizing beam splitter two or- 
thogonally polarized beams of frequency-nondegenerate 
squeezed vacuum. 

The generation process, depending on the experimen- 
tal conditions, is described by one of the Hamiltonians 



Hy± = ihg(a* H b v ± a v b ] H) 
H$± = ihg(a i H b^ H ± a v b v ) 



,t >,t 



h.c, 
h.c, 



(1) 



where g is the coupling constant proportional to the 
pump field, the PDC crystal length, and the second-order 
nonlinearity of the crystal. The down-converted photons 
carry linear polarization H (horizontal) and V (vertical) 
and are emitted in two frequency modes described by the 
creation operators a) and . The resulting states can be 
considered as macroscopic (multiphoton) generalizations 
of the two-photon polarization Bell states, 



^(oki^ia^&kHh.C.i 



(2) 



l^ic) = e r ^ b H±^)+h.c.| vac); 
where Y = J gdt is the parametric gain coefficient. 



The above equation allows one to determine their 

,(-) 



Schmidt decomposition. For the singlet state y&mac/, 
the decomposition is known For l^^ac), it 

similar form, so that both can be written as 
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l*^c>= E (±l) m \/K>^\n,m) a \m,n) b , (3) 



where X n = tanh 2 " Tj cosh 2 T and \n,m) a = \n) arl eg) 
\m) av denotes a two-mode Fock state with n photons po- 
larized horizontally and m photons polarized vertically in 
beam a (similarly for beam b). It is possible to factorize 
Eq. ^ further into two independent Schmidt decomposi- 
tions, one of them involving modes a#, by and the other 
one, modes ay,bu [28], 
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1 mac/ 



l*i >® 1*2), 



(4) 



E VKi\ n )a H \n) bv , 

n=0 

00 

£(±l) m V^|m)„ v |m) 6H . 
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one for \^>mac), but will be expressed in different polar- 
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it will be the right and left 
mlc) it will be the 



The Schmidt decompositions for the other two triplet 
states, |$mac)j can t> e easily written by recalling that 
they are obtained from j^mac) by rotating the polariza- 
tion. As a result, they will have the same form as the 

(+ 

rnc 

ization bases. For |<I> ma c 
(R,L) circular polarization and for |$ 
±45° linear polarization basis. 

The key difference between the two- and the four-mode 
BSV concerns entanglement. The two-mode BSV in- 
volves only photon-number entanglement. Namely, the 
photon number in the beam a with a fixed polarization 
(e.g. H) is unknown, but it is always equal to the photon 
number in the beam b having the orthogonal polariza- 
tion (V). This state is known to approximate the max- 
imally entangled EPR state in the high gain limit (30l |. 
In the coordinate representation in the limit Y — > 00 the 
electric field quadratures become completely uncertain 
and (5-correlated. Four-mode BSV is a product of two 
such states and thus, it simply provides two copies of 
it. However, it also contains polarization entanglement 
between beams a and b. This polarization entanglement 
is probed through measuring the photon-number corre- 
lations present in the two-mode BSVs (see Eq. J3])). It 
is called twin-beam multiphoton entanglement [24| and 



is most easy to notice if e.g. 
superposition 



14- 
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is re-written as a 
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*L»c> = 5- Y %/^TTtanh™ r\il>. 

mac cosh 2 r^ 1 



n 1 1 



(5) 
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where \tpn 



(-)\ 



particles. |\&iJ c 



is an analog of a singlet state of two spin-^ 

is invariant with respect to joint rota- 
tions of the polarization bases of both modes. Polariza- 
tion of each beam separately is undetermined, but the 
polarizations of the beams a and b are anti-correlated. 
The situation is similar for the triplet states. This ex- 
plains why these states can be considered as macroscopic 
(multiphoton) generalizations of the two-photon polar- 
ization Bell states. 



III. ENTANGLEMENT WITNESSES 

Entanglement witnesses are sufficient conditions for 
entanglement. Although sometimes being inconclusive, 
they are so far the only practical option for proving the 
entanglement of multidimensional systems. Conclusive 
(necessary and sufficient) conditions of entanglement are 
formulated only for special classes of quantum systems. 
For instance, the Peres-Horodecki criterion [8( provides 
a necessary and sufficient condition for the entanglement 
of two- and three-dimensional systems. Its continuous- 
variable counterpart, formulated by Simon is in the 
general case also only a sufficient condition, but becomes 
necessary for Gaussian states. In practice, continuous- 
variable entanglement is often witnessed using the Duan 
criterion containing variances of sum and difference 
quadratures for the subsystems of a bipartite system. 

In experiment, the output states l^mac), |$lfac) are 
given by a compound beam, comprising two independent 
frequency modes. Therefore, their Stokes operators Si 
are given by the sum of partial [3l| Stokes operators for 
modes a and b, 



S, = Sf + S° 



(G) 



with i = 0, 1, 2, 3, 6>Q = a^an + a) v av, Sf = a^an — 
Oyav, S 2 — a* H av + flyag, 6>f = i(ayan — aj^ay), and 
similarly for mode b. 

The condition involving variances of the Stokes opera- 
tors, 



Var(Si) + Var(5 2 ) + Var(5 3 ) > 2S , 



(7) 



holds true for any separable state of subsystems a and 
b 0, HH. This fact allows to formulate an entanglement 
witness operator 



W S = (S?+S b -(Sf + S?) ) 2 
+ (S^ + S b -(S^ + S b )) 2 
+ (S« + S b -(S« + S b )) 2 -2S . 



(8) 



As usual, negative mean value of the witness indicates 
entanglement. 

It is worth noting that the sign of (Ws) is invariant to 
the number of spatial and temporal modes because the 



overall state is a product ]T fe l v ^i ± ' ) ) °f states for differ- 



a-v.k, &H.k,bv,k and bn.k), so both the Stokes variances 
and the mean photon number of the whole beam con- 
tain additive contributions of separate modes. However, 
multimode separable states do not necessarily have this 
property. Therefore, the witness is valid for spatially 
and temporally multimode beams, under the assumption 
that separate modes are independent (the overall state 
is a product) I32|. 1331 . This property enabled its experi- 
mental testing [28] for the macroscopic BSV singlet state 

I^Lic). Indeed, the macroscopic polarization singlet Bell 
state has completely suppressed noise in all Stokes ob- 
servables S 1)2 ,3 E3, hence (Ws) = ~2{S ) < 0. 

At the same time, the witness will not be neg- 
ative for the three triplet states \^>m} c ), |$mac), since 
they have noise suppressed only in one Stokes observ- 
able [Hj]. However, all states in Eq. ^ can be trans- 
formed into each other by local polarization transforma- 
tions and thus, they contain the same amount of entan- 
glement. Based on Ws, we further derive the witnesses 
applicable to the triplet states. 

The l^irtac) state and the singlet state are linked by 
a local unitary rotation l^mac) = l a ®Ub\^mac), where 
Ub = e ,7rb ff hH , u\lAb = lb and l a ( fc ) is a unity operator 
acting on beam a (b) . This follows from the fact that the 
rotation l a <g>Ub transforms the Hamiltonian "H*- into the 
Hamiltonian %^ + . In experiment, this transformation 
is easily realized by means of a half-wave plate inserted 
into beam b with the optic axis oriented vertically or 
horizontally. Thus, the entanglement witness Wti for 
equals Wti = UbWsU\ , which yields 

Wti = (S? + S b - (S? + S b )) 2 
+ (S% - S b - (S% - S b )) 2 
+ (S%-S b -(S%-S b )) 2 -2S . 



* mac) 



(9) 



Of course, in theory the measurement of the witness Wti 
for the state l^^ic) is equivalent to the measurement of 
the witness Ws for the singlet state, (^~\Ws\^^) = 
(tf+|W T i|* + ). 
Entanglement witnesses for the other two macroscopic 

polarization triplet Bell states, |$^ a ' c ), can be easily ob- 
tained by recalling that the triplet states are transformed 
into each other via global rotations in the Stokes space. 

In particular, |<I>!mc) is obtained from |^^2 C ) by a tt/2 
rotation around the S 2 axis [33]. This rotation can be 
realized by a quarter-wave plate with the optic axis at 
angle 45° inserted into both beams a, b. The resulting 
witness can be obtained from ^ by changing variables 
Si — > S3 and 53 — > —Si, 



W T 2 = (5f - S b - (S? - S b )) 2 
+ (S% - S b (S% S b ) f 



S b )) 2 -2S . 



(10) 



Finally, |$L2c) 



is obtained from jvl/^c) by a 7r/2 rota- 



ent modes that are pairwise entangled (between modes tion around the S3 axis [33]. It can be realized by a tt/2 
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rotator (a quartz crystal or a Faraday cell). The witness 
has the form 

W T3 = (S? - S? - (S? - Si)) 2 



(SI 1 



S 



(S% 
(SI 



s b 2 )f 
s\)? 



2S . 



(11) 



Indeed, the witnesses (fMTTj) have negative mean values 
for the corresponding triplet states, because Var(£>f + 
Si) = Var(«S| - 5|) = Var(S£ -«S|) = for state l^ic), 
Var(5f-5f) = Var(5f-5|) = Var(5f +5|) = for state 



|$<+2 C ), and Var(«S? - «S X 6 ) 
5|) = for state 1*^). 

It is worth mentioning that the conditions for entan- 
glement given by witnesses (jSl fTTj) look very similar to 
the Duan condition as they contain the variances of sum 
and difference operators for subsystems A and B. At 
the same time, the relation between the Stokes operators 
from (|51 [Tl"j) and the quadrature operators from the Duan 
condition is only known for the case of light with a bright 
polarized coherent component and not for our case. 

In Fig. [1] we depicted an experimental setup where the 
witnesses Wt can be tested. It enables simultaneous 
measurement of various Stokes observables for beams a 
and b. The complete test will consist of the measure- 
ment of the variances for combinations of partial Stokes 
operators, Var(<Sf ± Sf ), i — 1,2,3, and the total pho- 
ton number (<Sq). The measurement should consist of 



Var(5 2 Q + S|) = Var(5 3 a 
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IV. ENTANGLEMENT MEASURES 

A. Schmidt number 

One of the well-known entanglement measures is the 
number of nonzero Schmidt coefficients \f\i in the 
Schmidt decomposition [2]. It is called the Schmidt 
rank or number. In case of a maximally entangled bi- 
partite system with symmetrical subsystems, defined in 
the Hilbert space H = C d (g> C d , all Schmidt coefficients 
have to be equal and the Schmidt number is K = d, 
where d is the dimensionality of the subsystem. In gen- 
eral, 1 < K < d and a state is separable if K = 1. To 
quantify entanglement in systems with infinite dimen- 
sional Hilbert space such as BSV, for which K = 00, 
another measure is more appropriate, which we will fur- 
ther call the effective Schmidt number. It is defined as 
follows i,S 



K^l/Tr(p 2 ) = 1/^2*1 



(12) 



where p is a density operator and ^ Aj = 1. This def- 
inition coincides with the original definition of K in the 
following way. For a separable state K = 1. For a maxi- 
mally entangled state with all Schmidt coefficients equal 
and d — > 00 (such a state does not exist because it is not 
normalizable) we would obtain K = d — >• 00. Otherwise, 
K is finite even if the number of Schmidt coefficients ^f\~i 
is infinite. 

The effective Schmidt number for l^mac) state is the 
product of the effective Schmidt numbers for subsys- 



tems I ^ 



1,2/ 



e: 



\2 
1 A,. 



which yields 



K = (1 + 2sinh T) 2 . Since all the macroscopic polar- 
ization Bell states ([2|) have the same form of the Schmidt 
decomposition, the effective Schmidt number for all of 
them is the same, K = (1 + 2A r o) 2 , where Nq = sinh 2 T is 
the photon population in each mode a#, ay ,bn,bv ■ For 
bright states Nq 3> 1 and K w 4JVq, hence the degree of 
polarization entanglement grows quadratically with the 
mean photon number. 



FIG. 1: (Color online) Experimental setup for testing the en- 
tanglement witnesses (f9l lll[) . In each beam, there is a Stokes 
measurement setup: PBS, polarizing beam splitter; HWP, 
half- wave plate; QWP, quarter- wave plate; D, detector. The 
signals from the two detectors in each arm are subtracted to 
obtain the Stokes observables. 

three series. In each series, one of the variances entering 
Eqs. (|9l llip is measured. This implies certain positions of 
the wave plates in the Stokes measurement setup. For S\ 
measurement, both wave plates should have their optic 
axes horizontal. For S2 measurement, the HWP should 
have the optic axis at 22.5° and the QWP, at 45°. Finally, 
for 53 measurement, the HWP should have the optic axis 
horizontal and the QWP, at 45°. The variances should 
be calculated by averaging over a large number of pulses. 



B. Negativity 

Other widely used and easily calculated measures of 
entanglement are the negativity and the logarithmic neg- 
ativity which gives an upper bound for distillable entan- 
glement j36[. For a bipartite quantum state p they are 
defined as Af(p) = Hp^Hi-l and Ej^(p) = log 2 \\p PT \\i, 
where PT denotes partial transposition with respect to 
one of the subsystems and ||»4||i = TrV A^A is the trace 
norm of an operator A. Particularly useful is the fact that 
the trace norm is factorizable, \\pi ® P2W1 — \\pi\\i\\P2\\i, 
and that the logarithmic negativity is additive, i.e., 
Ex(px ® P2) = E M { Pl ) + E M (p 2 ) d|. 

Since four-mode BSV is a product of two entangled bi- 
partite subsystems as in (U), it is sufficient to determine 
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the negativities of each of them, p\$,, separately. Per- 
forming partial transpositions with respect to the sub- 
systems by and respectively, we obtain ||pf T ||i = 
l|P2 T ||i = E^Lo V%] = e 2r - Thus, for the four-mode 
BSV, the negativity equals Af(p) — e 4r — 1. For high 
gain, N(p) ~ 16sinh 4 r = 16Nq. Again, the quadratic 
dependence of the degree of entanglement on the pop- 
ulation is observed. We notice that N(p) — AK. The 
logarithmic negativity takes the value Ejs/{p) = 4r/ln2. 
In comparison, for two- mode BSV it equals 2Y / In 2 (32J , 
which shows that BSV macroscopic Bell states contain 
(with respect to the logarithmic negativity) twice more 
entanglement than the usual two-mode BSV. This is un- 
derstandable: since four-mode BSV consists of two copies 
of entangled states, twice more entanglement can be dis- 
tilled from it than from a single copy. 



C. Fedorov ratio 

Both the Schmidt number and the negativity are not 
operational quantifiers of entanglement as they cannot be 
directly measured in experiment. For a bipartite system 
entangled in a continuous variable, an operational mea- 
sure was proposed [7|, called the Fedorov ratio. It is de- 
fined in the spirit of the entropy of entanglement with the 
advantage of being directly measurable in experiment. 
Further, we try to adapt this measure for characterizing 
BSV, entangled in the photon number, a variable that 
is discrete but can be viewed as pseudo-continuous when 
large numbers are involved. 

Consider a pure bipartite quantum system, entangled 
in a continuous variable with the state vector 

|tf) = J dv a du b F{v ai u b )\v a )\v b ). (13) 

The variable v a can be characterized by its marginal 
probability distribution P{v a ) — J dv b \F(v a , v b )\ 2 , with 
the standard deviation A^ a , and the conditional proba- 
bility distribution P{u a \v b ) = \F{v a: v b )\ 2 : evaluated for 
a certain value of v bl with the width 5v a . 

The effective Schmidt number K for the state jvP) is 
very well approximated by the ratio R Q defined as 



5v a 



(14) 



Of course, equivalently the variable v b may be involved in 
this definition, instead of v a . The parameter R, known 
as the Fedorov ratio, can be easily obtained in experi- 
ment and is hence an operational measure of entangle- 
ment. Note however, that this measure is only defined 
for pure states, so to be operational stricto sensu it has 
to be supplemented with an experimental proof of the 
purity of the global state. The possibility to generalize 
the Fedorov ratio to mixed states is an open question. 
For instance, it has been measured for the cases of wave 



vector [38| and frequency entanglement [3^|. For Gaus- 
sian states, the Fedorov ratio exactly coincides with the 
effective Schmidt number [Zoj . 

We adopt the definition (fT4")) for the photon-number 
variable and its probability distributions in the following 
way: 



R„ 



Sn a 



(15) 



where An a is the width of the marginal photon-number 
distribution in beam a, while 8n a is the width of the 
photon-number distribution in beam a under the condi- 
tion that a certain photon number n b has been measured 
in beam b. 

Similarly to the effective Schmidt number, the ratio R n 
for each of the macroscopic Bell states is a product of the 
ratios R n \ and R„2 of its two subsystems. For instance, 

for the state j^^ac), the states 2 ) m © have the same 
Fedorov ratios. They can be easily calculated by noticing 
that the marginal distribution P(n a ) is a geometric one, 



P(n a ) 



(tanhT 



\2n„ 



cosh^ r 



(16) 



while the conditional distribution is given by the Kro- 
necker delta, 

P(n a \n b ) = 8 na ,n b - (17) 
These distributions are schematically shown in Fig. [3J 
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: (Color online) Photon-number distributions for the 



Assuming that the width of the discrete P{n a \n b ) dis- 
tribution is unity, the ratios R n i,2 are given by the 
standard deviation of the geometric distribution (1161) . 

Rni,2 = V^sinh 2 r. Finally the Fedorov ratio for \^mac) 
equals 



2NJ 



(18) 



Thus, in case of high population, the operational measure 
R n differs from the effective Schmidt number K by only 
a constant factor 1/2. 
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The R n ratio for entangled BSV can be measured with 
the help of a setup shown in Fig. [TJ The orientation of 
the HWP and QWP should be such that proper polar- 
ization bases are chosen in the arms a, b. For instance, 
in the case of the singlet state \^ m ac}, the plates can be 
oriented in any way but similarly for arms a, b. The ra- 
tios R n i and R n 2 can then be measured independently 
using pairs of detectors D a \, D i and D a 2, D \. For each 
pair, after acquiring a certain (large) number of pulses, 
the photon-number distributions should be analyzed and 
the conditional and unconditional widths should be mea- 
sured. In practice, because the detectors do not distin- 
guish between close photon numbers, the photon-number 
distribution should be binned in the intervals of about 
200 photons [H|. 

All the above-considered measures of entanglement for 

the case of j^mic) states are plotted in Fig.[3Jas functions 
of the mean photon number. 

10 8 i 
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N 

FIG. 3: The degree of polarization entanglement in the log- 
arithmic scale for the four-mode entangled BSV state as a 
function of the average photon number No: the negativity N 
(solid line), the effective Schmidt number K (dashed line) and 
the Fedorov ratio R n (dotted line). 



D. Effective dimensionality of BSV Hilbert space 

Entangled BSV states of light are considered to be 
macroscopic generalizations of polarization singlet or 
triplet Bell states due to the symmetry reasons. Their 
polarization correlations are probed through photon- 
number measurements in polarization modes. Since the 
Hilbert space of these states is complex [4l[ and infinite, 
it is interesting to understand and quantify the amount 
of their accessible entanglement for example, by compar- 
ison with finite-dimensional systems, where the notion of 
a maximally entangled state is well understood. 

We propose a rough estimate for the dimensionality d 
of the Hilbert space of j^mac), depending on the gain T, 
based on the following argument. In two-photon experi- 
ments, small gain T < IO -3 is used to produce a super- 
position of the vacuum (the dominant component) and a 



biphoton. All higher-order contributions are largely sup- 
pressed. Truncation of the Fock states of order higher 
than one in Eq. ([3]), for a given gain, is justified if the 
normalization of the truncated state is preserved to a 
good approximation, given by a parameter e 

A 2 + 2A Ai = 1 - e. (19) 

For the two-photon case e « T 4 . However, it is known 
that in these experiments the four- and six-photon com- 
ponents are observed as well, which is manifested as, 
e.g., a decrease in the interference visibility for relatively 
bright PDC sources [42]. Thus, in our example, the value 
of e shows, for a given gain, how well the outcoming state 
from the PDC crystal can be approximated by a superpo- 
sition of a two-photon state and the vacuum. Of course, 
the smaller e is, the better the approximation gets. 

By analogy with the low-gain case, for any value of 
r we can locally restrict the state ([3]) to some finite- 
dimensional Hilbert space T-L — C da <g> C 6 , where d a and 
db denote dimensionalities of beams a and b. Then, 1 — e 
gives the probability to find the state in %. The natu- 
ral choice for the subspace C da is to keep in l^mic) only 
these components which have a limited number of pho- 
tons in beam a, a H ajj + et^ay < N max . Similarly for 
C db and beam b. This restriction implies the following 
normalization condition for the truncated state y&muc )'- 

N max N max -n 

J2 A « £ A '» = 1 - e - ( 2 °) 

n— m— 

The normalization is calculated over the sectors of the 
density matrix with fixed photon number, so that n+m < 
Nmax- It enables one to determine the dimensionality 
of C d °<», d = d a = d b = \{N max + l)(N max + 2), and 
the dependence of N max on the average population N . 
Using tanh 2 T = N Q /(N + 1) we turn Eq. fl^U) into 

e = (A) 1+Nma \N max + 2 - A (N max + 1)) (21) 

and obtain a linear relation between N max and No- It 
allows one to express the dimensionality for large popu- 
lation as d « - ^ Nq, where a is a function of e given 
by the equation e = e~ a (a+l). For example, if e = 10~ 12 
as is for the two-photon case, a s» 31. If e = 10 _2 (10 _1 ), 
we obtain a w 7(4). Please note the quadratic scaling of 
the dimensionality with Nq. 

The above estimations of the effective dimensionality 
of the Hilbert space for BSV are useful for reconstructing 
its density matrix, or the most significant part of it. It 
would provide almost full information about the joint 
photon-number distribution and could enable calculation 
of the entanglement measures based on the eigenvalues 
of the density operator obtained from the experimental 
data. 

Now, we investigate the amount of entanglement in the 
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FIG. 4: (Color online) K T /d as a function of e. The solid line 
is for No = 10, the dotted line No = 1, for large No (> 5) the 
curve does not change with No- 

truncated state 

n+m<N mam 

lC ( o t c ) ) = 7^s ^(±irv/AX:i«,m) a |m,n) 6 . (22) 

n,m=0 

Since the restriction is local, the amount of entanglement 
in \^mac) gives a lower bound on the overall entangle- 
ment in l^mac)- The effective Schmidt number for the 

truncated state K T = (1 - ef (E^i^** 

fulfills (j^) 2 K < K T < (1 - e)if . We look at how far is 

the restricted state y&mac ) from a maximally entangled 
state in % as a function of e, where H = C d <g) C d . This 
question may be understood in a more practical way as 
showing how good is our approach for creating maximally 
entangled state of dimension d. Since the maximally en- 
tangled states are defined with respect to dimensionality 
of their Hilbert space, we depicted K T /d as a function 
of e in Fig. 0J We notice that at high gain, a perfect 
maximally entangled state [K T j d — 1) is obtained only 

with e approaching one. This is because ySfmac) is a su- 
perposition of singlet states of different spin values and 
thus belonging to Hilbert spaces with different dimen- 
sionalities. Moreover, the genuine maximally entangled 
d-dimensional singlet \ip^ ) from Eq. J5]) has the smali- 

est weight in the superposition. Thus, the state \^>mac ) 
will certainly appear to be maximally entangled in de- 
tection, where it is projected on a fixed photon number 
subspace, but it is non-maximally entangled as a whole. 



An interesting fact is that the tradeoff between e and 
K T /d is independent of A^o above a certain value, which 
suggests that the "maximality" of entanglement in the 
discrete Hilbert space sense is unchanged when the gain 
increases. 

V. CONCLUSIONS 

In this paper we have discussed feasible entanglement 
witnesses and measures for bright squeezed vacuum in 
the form of macroscopic polarization Bell states, easy to 
generate in laboratories by means of parametric down 
conversion. Up to date, this is the only way to exper- 
imentally produce macroscopic entanglement (quantum 
correlations between two macroscopically populated sub- 
systems) for light. In general, this kind of entanglement 
is known to be very difficult to produce and to verify, 
since the more macroscopic it is, the more it is fragile 
to disturbances and the more it requires measurement 
techniques with very high precision. The last condition 
is difficult to fulfill with current-technology photodetec- 
tion. 

We derived the entanglement witnesses for the triplet 
BSV states and suggested their implementation based 
on the measurement of Stokes operators variances. This 
detection method is very convenient for intense fields be- 
cause it allows one to eliminate high intensity fluctuations 
present in other techniques, e.g. in direct intensity mea- 
surements. We also discussed the entanglement measure 
in the form of the effective Schmidt number, and found 
that it is very close to the Fedorov ratio, based on the 
measurement of conditional and unconditional photon- 
number probability distributions. We conclude that it 
can be considered to be an operational counterpart of 
the Schmidt number. Our results show that for all used 
measures, the BSV degree of entanglement scales as the 
mean photon number squared. 

We hope that the presented ideas will open the way 
for efficient entanglement verification for other infinite 
dimensional systems. 
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